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Abstract. We analyze the nonequilibrium dynamics of single inextensible semiflexible biopolymers as 
stretching forces are applied at the ends. Based on different (contradicting) heuristic arguments, vari- 
ous scaling laws have been proposed for the propagation speed of the backbone tension which is induced in 
response to stretching. Here, we employ a newly developed unified theory to systematically substantiate, 
restrict, and extend these approaches. Introducing the practically relevant scenario of a chain equilibrated 
under some prestretching force /pre that is suddenly exposed to a different external force /ext at the ends, 
we give a concise physical explanation of the underlying relaxation processes by means of an intuitive blob 
picture. We discuss the corresponding intermediate asymptotics, derive results for experimentally relevant 
observables, and support our conclusions by numerical solutions of the coarse-grained equations of motion 
for the tension. 

PACS. 61.41.-|-e Polymers, elastomers, and plastics - 87. 15. La Biological and medical physics. Mechanical 
properties - 87. 15. He Dynamics and conformational changes 



The multifaceted viscoelastic properties of complex cel- 
lular structures such as the cytoskeleton have been inten- 
sively studied on very different levels of complexity [1]. 
Concerning their nature as networks of semiflexible poly- 

■ mers, it has recently become apparent that prestressed 
[ networks provide a realistic model for the mechanical prop- 

■ erties of living cells [2]. Also, it has been found that the 
' dynamics of semiflexible polymers and their networks is 
, quite drastically altered by prestress [3-5] . On the level of 

■ single molecules, the scenario of a single semiflexible chain 
' prestretched with a force /pro that is suddenly exposed 
i to a different stretching force /ext at the ends may serve 
' to explain quite generically the influence of prestress on 
, single-polymer stretching and relaxation dynamics. Since 
' the flrst force-extension measurements on DNA, which cs- 
] tablished the wormlike chain (WLC) model as a very good 

theoretical description for semiflexible polymers [6], tech- 
nological advances towards signiflcantly improved time- 
and force-resolution [7, 8] have made measurements on 
such more involved experimental setups possible. Theo- 
retically, there exists a number of different approaches to 
the nonequilibrium dynamics of inextensible WLCs [9-12], 
most of them based on the concept of backbone tension, 
which is the force "holding the monomers together". Its 
key importance to the longitudinal dynamics has been rec- 
ognized some time ago [13], but so far its influence has 
mainly been treated on the level of heuristic scaling argu- 
ments [13-16]. While such reasoning based on the WLC 
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model has been used to analyze DNA relaxation experi- 
ments [17,18], a systematic theoretical description of prop- 
agation and relaxation of backbone tension has only re- 
cently been developed [19]. 

In the following, we consider an inextensible weakly- 
bending WLC in a simple solvent that is equilibrated un- 
der some tension /pre- At time t = 0, the chain is suddenly 
exposed to an external stretching force /oxt applied at its 
ends, the more interesting case being /oxt > /pre- Our ob- 
jective is to identify the dominant processes involved in the 
subsequent nonequilibrium dynamics, and, speciflcally, to 
determine how fast the contour stretches. The incxtensi- 
bility prevents any contour stretching in the bulk unless 
the end regions have been pulled apart, which is however 
limited by longitudinal Stokes friction. The conclusion is 
that initially the external stretching force penetrates the 
contour only within a boundary layer of size i\\{t), see 
fig. [TJ This comprises the relevant short-time effect and 
motivates the analysis presented in the remainder of the 
paper. After introducing model and equations of motion 
in sec. [11 we discuss in scc.[2]the stretching dynamics qual- 
itatively by means of an intuitive "blob" picture inspired 
by the analogon for flexible polymers [20]. Hereby we re- 
gard only the special case /oxt 3> /pro since it captures 
all prominent features of the stretching dynamics. We ob- 
tain crossover scaling laws for i\\{t) that are compared to 
literature results [13-16]. In sec.[3l we explicitly calculate 
tension proflles by solving the equations of motion for the 
asymptotic regimes of sec. [U and point out some incon- 
sistencies of previous approaches [13,16]. In the following 
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sec. m we present an algorithm to solve the equations of 
motion for the tension in intermediate and non-asymptotic 
regimes. Analytical and numerical results for experimen- 
tally relevant observables such as the change in the poly- 
mer's projected length are presented in the final sec. O 
where also some experimental implications are shortly dis- 
cussed. A more systematic derivation of the asymptotic 
equations of motion and their solutions, as well as details 
of the numerical approach are shifted to the appendices. 



1 Equations of motion 

Conformations of wormlikc chains arc described as contin- 
uous space curves r{s,t) of total contour length L. In the 
WLC-Hamiltonian [21], bending energy is proportional to 
the square of the local curvatureQ We set the proportion- 
ality constant (the bending rigidity) k = 1. Below, we will 
also set the friction coefficient (per length) C, to unity. In 
these units, tension is a length"^ and time a length^. The 
persistence length of a free WLC in equilibrium is then 
given by £p = (fcer)"^. We account for the local inex- 
tensibility |r'| = \dsr\ = 1 by introducing the backbone 
tension /(s, t) as a Lagrange multiplier function [9]. In the 
limit of small transverse displacements r± from a straight 
line, the corresponding weakly-bending Hamiltonian reads 



(1) 



Comparing contributions from bending T^//j^ and from 
tension fr\/t^ on the scaling level, we infer that on length 
scales smaller than the blob size 1^ = /"^''^ [24] (see 
fig- mil tfic conformation is dominated by bending forces, 
and only on larger scales perturbed by the tension contri- 
butions. In equilibrium, each blob carries a stored length 
of l^/ip [21], which is simply the thermal contraction com- 
pared to the straight conformation. The weakly-bending 
assumption requires ^ ip- This can be realized both 
for short stiff polymers (L ^ £p) and for strongly pre- 
stretched flexible WLCs (/pre > ^p^)- 

The overdamped motion of the transverse displace- 
ments is described by the equation dtr± = —SH/Sr±+^_^_, 
which captures the balance of viscous friction, bending 
and tension forces resulting from the energy [eq. ([T])], and 
stochastic noise 



dtr±_ = -rj 



(2) 



Introduced as Lagrange multiplier function, the tension 
/(s, t) has to be computed through the inextcnsibility con- 
straint, which makes eq. ^ highly nonlinear. In spite of 
that, it can be simplified quite significantly in the weakly- 
bending limit [19]: spatial variations in the tension / can 
be neglected. It turns thus into a linear equation that can 
be decomposed into normal modes. 



dtr± = -q^r±_ - q\fr±_ + ^ 



(3) 



^ Note that this harmonic dependence has lately been under 
some debate concerning DNA bending angle distributions on 
short length scales [22,23]. 



and solved using the response function 

X±{q;t,t') =e-«'[«'(*-*''+/t'^'^*]6>(t-t')- (4) 

Here, 0{t) is the Heaviside function. This allows to eval- 
uate a characteristic quantity of high relevance for our 
approach: the stored length density Q=\ {''''T) that mea- 
sures the average amount of contour length stored in ther- 
mal undulations: 



E U(9)xi('z;i,0) + |^ / dt'xl{q-Xt') 



(5) 



by evaluating the noise correlation fc(0^ j. g(i')) = 
4/ {L£p)6k^qS{t — t'). The initial mode spectrum of a poly- 
mer equilibrated under the force /pro follows via equipar- 
tition from the Hamiltonian [cq. 



Qo{q) 



LW + /pr 



(6) 



However, as discussed in the introduction and illus- 
trated in fig.[l] neither the tension / nor the stored length 
density g are in fact spatially constant. Compared to the 
"fast" transverse undulations though, their spatial depen- 
dence is rather slow and can thus be re-introduced within 
an adiabatic approximation. This convenient practice is 
rigorously justified by means of a multiple scale analy- 
sis in rcf. [25]. Its main result is that the curvature in 
the spatially slowly varying tension profile is given by 
changes in the stored length density g{s,t), which inher- 
its its slow arclength dependence adiabatically from the 
tension f{s,t) via eq. ([4]): 

d'j{s,t)^-dtg{s,t). (7) 

Integrating cq. ((T]) over time yields 

d^,F{s,t)^~[gis,t)^g{s,0)], (8) 

where we have introduced the integrated tension 



F{s,t)= f dt'f{s,t') 
Jo 



(9) 



Using eq. ([5]) with a spatial dependence of g{s, t) given 
through /(s,i) in the exponent of x±iq'jtit')i E^nd taking 
the continuum limit L —> oo, renders eq. ([8]) in the explicit 
form [25] 



' dq f 1 - e-2<z'[9't+^'(s.t)] 







pre 



-2q^J dt'e-29'[9'(*-*')+^(^^*)-^(^.*')]|. (10) 

Initial and boundary conditions for our setup are 
F(s, i < 0) = /prct and ^^(0, i) = F(L, i) = /cxtt. (H) 
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Fig. 1. Pulling a prestretched filament. Thermal fluctuations 
in the contour r(s) are straightened first in a boundary layer 
of width £\\{t). At any time t, blobs of size fprl^^ and /"^''^(s) 
are associated with the prestretching force /pro and the actual 
local backbone tension f{s), respectively. 



For a detailed derivation of eq. (|10|1 we refer the in- 
terested reader to ref. [25]. Here we only comment on a 
few important points. First, a proper decomposition of 
eq. ([2) into its eigcnmodes depends on the choice of bound- 
ary conditions and yields, in general, eigenvalues different 
from the simple q^[q^ + f] used here. Second, we are going 
to extend eq. (fTO|) . which applies to the continuum limit 
L — !■ cxD, to "real" polymers of finite length. Both issues arc 
resolved by recognizing that on short times the relevant 
modes have wavelength much shorter than L, that long 
wavelengths are initially suppressed by prestretching with 
/pro, and that therefore the influence of the boundary con- 
ditions and the discreteness of the mode spectrum is small. 
Below, we will discuss the relevance of these effects for 
pertinent observables. Finally, the underlying scale sep- 
aration between small-scale "fast" transverse fluctuations 
and large-scale "slow" tension dynamics can be readily ra- 
tionalized by observing that the tension varies only over 
distances on the order of ^|| , which is at any time much 
larger than the correlation length for transverse fluctua- 
tions [19]: 




for t < /"^ 
for t > /-2 



(12) 



Eq. (fT2|) is easily understood as the expression of the 
dynamic force balance between transverse friction forces 
r±/t and bending r^/f^ or tension terms fr±/i'^, respec- 
tively, i.e., by analyzing eq. ([2]) on a scaling level. It dis- 
plays the crossover between "free" (bending-dominated) 
relaxation within blobs and "forced" (tension-driven) re- 
laxation on larger scales. 



2 Blob picture of stretching dynamics 

A qualitative yet rather thorough understanding of the 
noncquilibrium stretching dynamics can be gained already 
from the special case /ext ^ /pre, i-e., for the scenario of 
a sudden strong increase in stretching force. It is useful to 
interpret the correlation length £± (t) of eq. (fT^ as an equi- 
libration length for transverse fluctuations. At any time t. 



segments of length £j_ arc in equilibrium with their sur- 
roundings [15, 19]. Upon applying the external force /ext, 
the contour starts to stretch within a boundary layer of 
size Decomposing the latter into £\\/£± segments of 
length £±, the extension S of each segment in response 
to the local tension f{s) can be estimated within equilib- 
rium theory. At this point wc use the blob picture: the 
prestretching force /pre and the local tension f{s), which 
builds up after the external force is applied, induce blobs 
at different length scales. Associated with the weaker pre- 

— 1/2 

stretching force are large blobs of constant size /pre , and 
the stronger local tension corresponds to small blobs of 
varying size /~^/^(s), respectively (see fig.H]). The exten- 
sion S depends on "how many" blobs are contained in a 
segment (or vice versa). Three cases can be distinguished 

for /ext > /p rc- 

^ ^-L ^ f~'^^^{s): a large number of segments £j_ are 
in either type of blob; hence, bending forces dominate 
and the extension follows from linear response [26]: 

5 ~ i'if{s)/ep. 

_ j-i/2^g-j ^ ^ /prc^^: segments £j_ near the bound- 
ary arc larger than the small blobs corresponding to 
the local tension. Hence, they get almost completely 
stretched, and essentially all initially stored length is 
pulled out. Being still smaller than the large blobs of 

size /pr^^, these segments correspond to short stiff ini- 
tially unstretched chains, thus their extension is (5 ~ 
£l/£^. 

— 1/2 

^ /pre ^ £±' now the segments i± are larger than any 
of the blobs and release the stored length of a taut 
string of blobs. In a segment, there are ^xZ/prc^^ large 

— 1/2 

blobs of size /pre each with stored length l/(^p/prc), 
hence S ~ £±/i£pfp(c)- 

The above results for the extension 6 and eq. (fT2|) for 
the length £± of the segments can now be used to estimate 
the total stretching (^|| /£±)6 of the chain's boundary layer. 
On the scaling level, we can set /(s) = /ext for the relevant 
segments near the boundary. The size £|| of the boundary 
layer is the central quantity still to be determined. Ac- 
cording to its definition, it is obtained by requiring that 
the total longitudinal friction on the order of£^^{£^^/£±)6/t 
equals the driving force /ext, i-e., £\\it) scales as 



£\/V/\ 



for t < fj 



(13a) 



^;,/'(/ext ^)'/^ for /,-? « t « (/ext/pre)-' (13b) 

4^'/p/e'(/ext0'/', for (/ext/pre)-' « i. (13c) 



The first case [eq. (|13ap ] has already been derived by 
Everaers et al. [15] (EJAM) based on a similar argument. 
This universal initial regime also shows up for other force 
protocols [19,27], including the case where /ext < /pre- 
The second scaling law [eq. ()13bp ] has first been proposed 
by Seifert, Wintz, and Nelson [13] (SWN). This interme- 
diate regime emerges if the large and small blobs have 
very different size, i.e., either for vanishing prestretching 
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Fig. 2. Scaling of the boundary layer size i\\{t) (log-log 
scale) for different force ratios /prc//cxt = 1/2 (long-dashed), 
/pre//ext = 10"'' (short-dashed) , and /pro ~ (solid). £\\ is ex- 
tracted from numerical solutions of eq. (|10|l . The asymptotic 
scaling laws of eq. (|13p are indicated by dark lines. 




Fig. 3. Schematic tension profiles for /oxt > /pre. Left: propa- 
gation regime t ^ t^. Near the ends, the tension falls off over 
the distance ^ L and is constant in the bulk. Right: re- 
laxation regime t ^ t\. The tension relaxes towards its final 
equilibrium value / = /ext. For times t ^ the difference 
between / and /ext is small. 

i^/^-regime of eq. (jl3bp appears for very strong force scale 
separation only. We proceed with a more detailed discus- 
sion based on tension profiles to elucidate the similarities 
and differences between our results and previous work. 



force or for strong force scale separation /cxt ^ /pre. How- 
ever, only near the boundary the different blobs differ in 
size. Accordingly, the contour stretches significantly only 
where the actual tension is already strongly different from 
its bulk value /pro- Thus, as pointed out in ref. [14], ten- 
sion propagates faster than stretching is achieved. Finally, 
for long times t ^ (/oxt/prc)~^, we get a result similar 
to the one by Brochard, Buguin, and de Gennes (BBG) 
derived in ref. [16], where they assumed the tension to 
be locally equilibrated. Note, however, that we obtain a 
more detailed force dependence in eq. (|13cp . In sec. |3] be- 
low, this discrepancy is traced back to a subtlety of the 
limit /pro 0. Both results agree only for /ext « /pre- 

In the latter case, i.e., if the stretching force is changed 
only by a small amount Af = /cxt - /pro < /ext, the re- 
sults for the segment extension S are reduced by a fac- 
tor Zi///oxt- Because the different blobs have then es- 
sentially the same size, the force balance argument lead- 
ing to eq. (here with the driving force Af) yields 

only two cases: £|| ~ il^'^t^^^ [15] for t < /"? and i\\ ~ 

^y^/ext*^^^^ [16] for t ^ /cxt- There is no intermediate 
t^/^-regime. 

Summarizing the preceding discussion, we found inter- 
mediate asymptotic scaling laws for the boundary layer 
size £|| (<) by estimating the stretching of segments £j_ 
within a blob picture and balancing the resulting longitu- 
dinal friction with the driving force. As a rule of thumb, 
the small blobs corresponding to the stronger force /(s) 
decide whether a segment gets stretched "just a little" or 
"much" , and the large blobs corresponding to the weaker 
force /pre determine how much "much" actually is. To 
illustrate the power laws given in eq. p^ . numerical re- 
sults for €|| (t) are shown in fig. [2l where (t) has been 
extracted from tension profiles computed by numerically 
solving eq. PU)) . see sec. SI The intermediate asymptotic 
scaling is clearly visible for small and large times, and the 



3 Tension profiles 

Instead of repeating the rigorous analysis of eq. pI7|) pre- 
sented in ref. [27] and extending it to deal with two force 
scales, we will analyze three different asymptotes of eq. PH)) 
motivated by physical considerations that lead to the quite 
different approximations employed previously [13,15,16]. 
As in sec. [21 we are primarily interested in the case /oxt 3> 
/pre and discuss other possibilities only incidentally. A 
more systematic and complete derivation is presented in 
appendix lA.ll Having motivated these approximations and 
derived their asymptotic validity, we will find that the par- 
tial integro-differential equation (PIDE) [eq. pU)) ] reduces 
to simpler differential equations that can again be solved 
in two regimes (see fig. [3]): 

1. In the propagation regime t <C t| with t| defined via 

^11 (^1) ~ ^' tension profile falls off over the (in- 
creasing) distance £\\{t). Because then £^^{t) ^ L, we 
may neglect the presence of a second end, set L ^ oo 
and work with a so-called "semi-infinite" contour. For 
this regime we will show that, using SWN's and BBG's 
approaches, the resulting differential equations do not 
always lead to tension profiles with the required fea- 
tures (compare fig. namely that F{s) decays from 
/oxt t (at the boundary) to /pro t (in the bulk) within 
a region of size ^y , and that F{s) = /pre t is flat in the 
bulk. 

2. In the relaxation regime t ^ the influence of the 
second end cannot be neglected. However, it turns out 
that for late times < ^ t*, the resulting tension pro- 
file can be obtained by expanding eq. pU)) about the 
final equilibrium profile: F(s,t) — fcxtt + 5F{s,t). 
To lowest order, the correction SF is then a simple 
parabola with constant curvature given by the right 
hand side of eq. pU]) evaluated with F = /cxt t. The 
time ti, is the crossover time to a regime t ^ t^, where 
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6F <C fcxt t is in fact small, and can therefore be de- 
fined via 6F{ti,) ~ fcxtt*- With one exception, we will 
find that ~ . Some of the relaxation regimes have 
already been discussed in refs. [19,27] for related force 
protocols. We shift their analysis to appendix lA. 21 and 
consider in the following only generically new results. 

The scaling of the crossover times <| and ti, between the 
respective propagation and relaxation regimes will be de- 
rived along the way and is summarized in table [T] below. 
It is different for the three approximations analyzed in the 
remainder of this section and depends on the magnitude 
of /pre and /ext relative to the critical force /c = ip/L"^- 
Similar to the well-known Euler buckling force /* oc L^^ 
that represents a critical threshold above which forces sig- 
nificantly disturb the transverse dynamics [3,19], the force 
fc can be seen as threshold force for the longitudinal dy- 
namics. 



reads 



Fis,z) 



Z\/ cosh[(L - 2g)/2A] /pre , . 
cosh(L/2A) z2 ' ^ ' 



where Af = /oxt — /pre is the force difference. 

If L ^ A, the tension profile eq. pT)) varies only close 
to the boundaries, as it is characteristic for the propaga- 
tion regime. Near s = 0, it simplifies to 



F{s,z) 



-s/\ 



/pr 



Z^ Z^ 

which can be backtransformed [27] to 



Fis,t)^Aft<p{s/£^\{t)) + fp,,t, 



(18) 



(19) 



where 0(^) w exp[— 2 '^/^^/F (15/8)] is a scaling function 
that depends only on the ratio ^ = s/^||(i) with a bound- 
ary layer size i\\{t) = ip^'^t^^^ that scales as predicted in 
eq. (|13ap . The requirement L ^ X translates into t <C tf^ 



with t 



3.1 Linear approximation 

The two limits discussed above arise quite naturally within 
the linear approximation to eq. pop . which is similar to 
the approach by EJAM [15]. It consists in treating the 
prestretching force /pre as well as the tension F as small 
perturbations with respect to the dynamically relevant 
bending contributions. From eq. (fT2|l and its interpreta- 
tion we infer that this approximation is justified in the 
limit t <|; f~^- Because the tension / varies between /ext 
at the boundary and /pre in the bulk, this certainly holds 

if t < min(/-?, /-2). In this case, we approximate the re- 3.2 Taut-string approximation 



The relaxation regime corresponds to the complemen- 
tary limit L ^ A, or equivalently i ^ Here we find, as 
anticipated, that the tension F{s, t) differs for t ^ ti, from 
the flat profile /ext t only by a subdominant term dF ~ 
/exti^^*i^/^p with a (trivial) parabolic s-dependence. By 
evaluating SF{t^) = f^^tU we find t^ ~ L^/f^ ^ ^11 
appendix lA.2l for details. 



t'l. See 



spouse function [eq. (g])] by x±(g;t,i') « e'^"? *{l-2q'^F) 
and the initial mode spectrum [eq. ^] by {q^ + fpic)^^ ~ 
q~^{l — fpre/q"^), and obtain eq. ([T0|) to linear order in F 
and /pre as: 



«^<"''"=y„ lit 



dq 



/pre 



1 - e 



+ 2F(s,t)-V / dt'F(s,t')e-29'('-*') 



In this approximation introduced by SWN [13], bend- 
ing and thermal forces are neglected after the prepara- 
tion of an initial equilibrium configuration. As discussed 
below eq. (fT^ . bending contributions are locally negligi- 
ble against tension in the long-time limit t ^ /^^. The 
transverse displacement modes effectively obey dtr±^ = 
^<ffi'±^ see eq. The stored length density relaxes as 
Q = Y.q Qo{q) exp[-2q'^F], and eq. 0} reads [27]: 



(14) 



Using the Laplace transform F{s,z) = C{F{s,t)}, this 
reads 



d'^F 



dq 1 



TI'^P + /pr 



(20) 



dlF{s,z)^ 



2 /pre 



Depending on the magnitude of the product fp^cF (which 
is dimensionless in our units), this has two asymptotes: 



j-+F(s,z) 



2z 



(15) 



z(z + 2g4) ' ^ z + 2q'^ _ 

which, after performing the g-integral, reduces to: 

X^d^J^F^^. (16) 



diF 




2F/7t£1, 

l/(2/pVe^p), 



ff fprcF < 1, (21a) 
ff fp,,F » 1. (21b) 



Here, A(z) = 2^^^£p^^z is a dynamic length scale de- 
noting the size of spatial variations in F{s, z). It is there- 
fore directly related to ^|[(t). With the boundary condi- 
tions F(0, z) = F{L, z) — /ext/z^, the solution to eq. (fT6] 



These will be discussed only for the propagation regime. 

Estimating F ~ /exti near the boundary, the first 
asymptotics is realized for intermediate times t 
(/pre/ext)~^- The differential equation (|21ap is solved [27] 
by the scaling ansatz F{s,t) = fcxtt (t>{s/£^^{t)), with the 

boundary layer size £\\{t) = £p^^ (fcxtt)^^^ as in eq. (jl3bp . 



which implies t j^ 



1/2^ f ,^^l/4 

L'*/(£p/ext)- This scaling ansatz gives 



6 



B. Obermayer et al.: Stretching dynamics of semiflexible polymers 



a tension that decays within a length £\\(t) as expected, 
but does it meet our expectations in the bulk? If /pro = 0, 
corresponding to the "pulling" -scenario of ref. [19], the 
scaling form F = /oxt^ with ^ = s/£\\{t) can smoothly 
be extended to the bulk tension, which obeys F = d^F = 
[27]: 



1 



(727r)i/4 



e (727r) 



1/4 _ 



(22) 



If /pre is finite, however, the bulk tension has nonzero mag- 
nitude and cq. (j2ip would yield nonzero curvature for the 



expectedly flat bulk profile. Hence, for finite /pro the taut- 
string approximation is not valid along the whole contour, 
but only in a finite region near the boundary. In this case 
it does not represent an intermediate asymptotic regime of 
eq. (Uni). 

The second asymptotics (j21bp . realized for late times 
^ ^ (/pre/ext)~^, implies zero curvature for the tension 
f = dtF everywhere. The relaxation changes its character: 
most of the initially excited modes have relaxed (the long- 
wavelength contributions have been "cut off" by the pre- 
stretching force from the beginning). This corresponds to 
an almost completely stretched contour under linearly de- 
creasing tension in the boundary layer. However, eq. (j21bp 
is not sufficient to describe the specific shape of the asso- 
ciated tension profiles which have to be constant in the 
bulk. One needs to include thermal noise0 



3.3 Quasi-static approximation 

The preceding paragraph revealed a problem with the 
taut-string approximation for f^t 3> 1 and fpreF 3> 1: 
the right hand side of cq. (j21b[) becomes time-independent 
and the deterministic relaxation saturates, so that the sub- 
sequent dynamics is purely of stochastic origin. This al- 
lows to simplify it as a quasi-equilibrium process. Follow- 
ing the approach of BBG [16], one can assume the fila- 
ment to be equilibrated under the local tension /(s, t) = 
dtF{s,t). From the continuum approximation to eq. 
the stored length density is g = J ^{q'^ + dtF)~^ = 

(9tF)-i/2/(2£p). According to cq. ©, this adds a small 
but relevant contribution to eq. (j21bp [27]: 



dlF 



1 

2L 



-1/2 



By taking a time derivative we get 

dtf 



4£p/3/2- 



(23) 



(24) 



We discuss the quasi-static approximation both for the 
propagation and the relaxation regime. 



SWN's "noiseless" simulations [13] for a prestretched chain 
yield the correct scaling £j| (t) oc t^^'^ and linear tension pro- 
files, but (as their scaling argument indicates) a straight string 
drawn through a viscous solvent gives a linear profile, irrespec- 
tive of thermal noise. 



1000 
4 




Fig. 4. Numerical solutions of eq. (|25|l for the scaling func- 
tion iy3(^) of the tension profile f{s,t). (a) Plot of the slope 
d^ip{0) at the boundary (solid line) vs. force ratio c — 
/prc//oxt (log-scale) with analytical asymptotes for c — > and 
c ^ oo, respectively (dashed lines), (b) The scaling func- 
tions (c = 2, 4, 10, 20, 100, 1000 from top to bottom) asymp- 
totically collapse for c — > oo if the abscissa is rescaled by 
c~^^^ and the ordinate by c~^. (c) The scaling functions 
(c — .5, .25, .1, .05, .01, .001 from top to bottom) asymptoti- 
cally collapse onto a piecewise linear profile for c ^ if the 
abscissa is rescaled by c~^/*. 



Propagation regime. Here we insert the scaling ansatz 
f{s,t) = fcxtfiO with ^ = s/£jj(t) and the tentative scal- 
ing q(t) = el^^f^l^t^^^ as proposed by BBG. This leads 
to an ordinary differential equation 



(25) 



Boundary conditions arc ip{0) ~ 1 and (p(^ —s- oo) = c 
where c = /pro//cxt is the force ratio. This is the same 
equation as eq. (18) of ref. [16], safe for a factor j due to 
a slightly different definition of the scaling variable ^. We 
have solved it numerically by a shooting method. Starting 
with ip{0) = 1, it is integrated forward, adjusting the slope 
d^(p{0) at the left boundary in order to fulfill the second 
boundary condition. The dependence of the slope d^ip(0) 
on the force ratio c is shown in fig. |D(a), together with 
the asymptotes for c — *■ oo and c — *■ 0, respectively. These 
shall now be analyzed in more detail. 

The case c > 1 corresponds to a stretching force /oxt < 
/pre that is decreased at t = 0, i.e., a "sudden release"- 
scenario as in refs. [16,19,27]. Since /pre is the relevant 
force scale, the scaling ansatz is more appropriately writ- 
ten as f{s,t) = c/extv(c~'^^^0- ^ow we can safely take 
the limit c — > oo (/ext ~^ 0) describing the case where 
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the force is switched off completelyH because the artificial 
c^'^^-divergence in d^ip{0) has been removed. The resulting 

boundary layer size scales like i\\{t) ~ ip^'^ fplot^^'^ as in 
ref. [16] and the corresponding scaling functions, shown in 
fig. [UJb) , smoothly converge to the solutions depicted in 
refs. [16,27]. 

The opposite limit c — > is more subtle. Although 
it is not obvious from eq. (j25p . we may not simply set 
/pre = as in ref. [16], since then the condition fprcF ^ 1 
in eq. (j21bp cannot be met. Moreover, considering the 
limit c — *■ with /pro > fixed instead, our numerical 
results (fig. mja)) indicate that the initial slope d^ip{0) di- 
verges like c~^^^. This suggests the scaling variable rj = 
= s/mt) with the new boundary layer scaling 

^ll(t) = c^/^qit) = 4^Vp/o'(/oxtt)^/^ as anticipated in 
eq. (|13cp . We can further rationalize this asymptotic scal- 
ing by inserting the improved scaling ansatz F(s,t) = 
/exti(/>(77) into eq. which gives d^^(j){r]) = \ + 0[c^l'^). 
In the limit c — > 0, the scaling function is parabolic in 
the boundary region r\ £ [0,2], and the scaling function 
= — ^rjdtjcj) =1 — ^77 for the tension / becomes linear. 
Numerical solutions to eq. (|25p in terms of 77 = ^c~^/* in 
fact tend towards a piecewise linear function with a kink 
fixed at 77 = 2 (see fig. Hljc)), hence the kink is moving 
towards vanishing ^ = as c — > 0. This reveals that 

the limit c —s- cannot properly be taken in terms of the 
original scaling variable ^, because eq. ([23| and therefore 
eq. ([25]) become invalid if c = /pre /./oxt = 0. 



Relaxation regime. In the relaxation regime t ^ ti, we 
expect the tension to deviate only slighly from the flat 
equilibrium profile F = /oxt t, see fig. El We therefore write 
F{s, t) = fc^tt + SF{s, t) and expand eq. ([23]) . which gives 
the time-independent correction: 



6F 



1 



-1/2 



pre 



,-1/2 

J cxt 



s(s-L). 



(26) 



This result has two significant implications. First of all, in 
order to estimate the time beyond which 5F becomes 
small compared to /oxt t^ we evaluate 5F{ti,) ~ /oxt t* and 
find that not always ti, ~ i|, see the last two rows of 
table [T] We find for /pro ^ /oxt and /pre ^ fc another 

intermediate regime <C i <C of homogeneous ten- 
sion relaxation, where the tension / has a more compli- 
cated than just parabolic spatial dependence and the sim- 
ple ansatz F{s,t) = fcxtt + SF{s,t) for the relaxation 
regimes does not work. This has been attributed [19] to 
the fact that longitudinal friction may play an essential 

role also for times i 3> t| beyond the propagation regime 
and significantly slow down the relaxation. Safe for small 



If for strongly prestretched "flexible WLC" with £p <^ L 
and /pro ^ the external force is decreased too much 
(/ext ^ ^p'^)i the weakly-bending approximation eventually 
breaks down near the ends as the tension decreases. This re- 
stricts the applicability of our analysis to the dynamics of those 
observables that are insensitive to boundary contributions [27] . 





4 


t* 


/pre <IC /ext ^ fc 




/cxt ^ /pre ^ fc 


/pre fc /cxt 




fc /pre /cxt 


-^^ / (^P ./pre /ext ) 


/ext fc /pre 


L /(-^p/pre ) 




fc /ext <^ /pre 





Table 1. Crossover times and depending on the magni- 
tude of /pro and /cxt relative to the critical force fc — ^p/L*. 



correction terms (see appendix IA.2P , our results for the 
homogeneous relaxation regime are identical to the "sud- 
den release" -case /ext = discussed in ref. [27]. 

Secondly, a time-independent correction SF of the in- 
tegrated tension yields an actual tension / = dtF that 
(to leading order) does not deviate at all from the flat 
equilibrium profile / = /ext- The relaxation dynamics sat- 
urates prematurely to equilibrium, i.e., at times t ti, 
long before the longest bending mode has relaxed (see 
also eq. ((37)) below). This behavior arises within a sce- 
nario where a strong stretching force /pre ^ fc is changed 
to /ext ^ fc but remains strong. Hence, there is no need 
to destroy (create) long-wavelength contributions in the 
mode spectrum of the stored length (see eq. ([6])) as in 
the "pulling" ( "release" )-scenario of ref. [19] . By taking the 
next-to-leading order term in the quasi-static approxima- 

wc find an improved cor- 



tion [eq. (cf. appendix IAT2|) 

rcction term that captures the slow diffusive relaxation of 
the remaining long- wavelength modes: 



5F{s,t) = 



1 



f-1/2 _ 

J pre J cxt 



,{s - L) 



f pre / e: 



4£p V2^(/extt)l/2 



s(s-i). 



Having discussed some of the asymptotic solutions to 
eq. (fTO]) resulting from the hnear, the taut-string, and the 
quasi-static approximation, respectively, we now present 
a numerical approach to solve eq. ([TI 



4 Numerical approach 

While the intermediate asymptotes of eq. PU)) serve to 
expose the relevant physics involved in the relaxation pro- 
cess, experiments on the biopolymers we have in mind 
are usually performed in intermediate, non-asymptotic pa- 
rameter ranges. Therefore we devised a numerical scheme 
to solve the nonlinear PIDE [eq. (fTO|) ] under very general 
initial and boundary conditions, allowing cq. (jlOp to be 
applied to far more scenarios than just the ones analyzed 
in this work. 

Numerical solutions of cq. (jlOp are written in terms 
of dimcnsionless variables. Here, we use a characteristic 
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length scale Sc specified when necessary, and the time and 
force scales tf = and /ext, respectively. The integrated 
tension is written as 



(27) 



and the actual tension profile is then extracted via 

/(s,t) - /ext(/^(a,T) = f,^tdr[T${a,T)]. (28) 

Inserting the scaling form [eq. ([27]) ] into eq. (fTO|) and 
rescaling q = /oxt and t' = zt, we are left with 



'dq 



1 



T{q^ + c) 

-2f I d2e"^'''^t'^'^^"^'+*('^'^)~'^*^'^'^'^)l 
Jo 

The remaining parameters are 

= £p/ife!t4) and c = /prc//cxt. 



(29) 



(30) 



For the "semi-infinite" polymers with L —> oo we use = 
1 by choosing Sc = £p^'^ f^J^.^^ , while in the finite case the 
appropriate choice is Sc = i. 

Our strategy is as follows. Given the initial condition 
(p{a,0) = (po{a) and appropriate linear boundary condi- 
tions for the respective scenario, we introduce a discretized 
time coordinate t„ for n ^ 0, . . . N . At each time step t„ 
we obtain a two-point boundary value problem: 



dl^ia,Tn) ^ G[^ia,Tk)k<n]. 



(31a) 



The nonlinear term G depends also on "earlier" solutions 
^(ct, T;;) with k < n. Wc allow arbitrary linear boundary 
conditions incorporated via 6 coefficients ajk'- 

aoo^(0, T„) -t- aoida^iO, r„) = ao2 (31b) 
aiQ<P{aM,Tn) + anda<l'i<JM,Tn) ai2. (31c) 

These boundary value problems [eq. (|3T|) ] are solved by 
transforming to a system of nonlinear equations using the 
discretized arclength coordinate am for m = 0, . . . M , and 
discrete representations for the differential operators da 
and For details see Appendix [BJ 

Given numerical solutions to eq. pU]) , it is straightfor- 
ward to obtain results for pertinent observables. 



5 Results for the change in projected length 

The observable of highest interest in experiments is cer- 
tainly the change in end-to-end distance. Since the sudden 
change in stretching force (from /pre to /ext) induces the 
creation or destruction of stored length, and since the to- 
tal contour length is conserved, we can simply evaluate the 



change in end-to-end distance (approximated by its pro- 
jection onto the longitudinal reference axis) by integrating 
up the total difference in stored length: 



Z\i?ll(t) = - / ds[g{s,t)- g{s,0)] 
Jo 

:)2 7 



<lsdiF{s,t) 



(32) 







F'{L,t)-F'{Q,t). 



Here we have used eq. © to relate the change in stored 
length to the curvature of the integrated tension F. The 
notation AR\\{t) serves as a remainder of the fact that 
Z\i?ll {t) includes only bulk contributions. The "real" change 
in projected end-to-end distance (/4_R||(t)) measured in 
experiments also contains end contributions that depend 
on the choice of boundary conditions which we have ne- 
glected in our approximate mode decomposition of eq. ^ . 
However, a careful analysis [27] shows that for hinged or 
clamped ends and a sufficiently large prcstretching force 
/pro ^ these contributions vanish. For the "semi- 
infinite" polymers discussed above we neglected the pres- 
ence of a second end; here we can exploit the symmetry 
of the setup and use Z\i?||(t) = -2F'{0,t). 



5.1 Analytical results 

Given the two forces /oxt and /proj and the additional force 
scale /o, the sequence of asymptotic regimes realized for a 
specific choice of parameters depends on their respective 
ratios and is in general quite complicated, see also tabled] 
Relaxing the condition /ext > /pre of sec. [2l we chose to 
first illustrate schematically the general dependence on 
the force ratio /pre//ext since the systematic investigation 
of its influence is the major novelty of this work. With 
our numerical results, we will then demonstrate the de- 
pendence on the polymer length (more specifically, on the 
parameter of eq. (|30p) for fixed force ratio. 

Using eq. (j32p , we summarized the asymptotic growth 
laws for (t) resulting from the tension profiles com- 
puted in sec. [3] in two phase diagrams. We illustrate the 
growth laws for AR\\ (t) as boxed formulas in their re- 
spective asymptotic regime depending on the force ratio 

/pro//oxt in fig. [5] for the propagation regime (i ^ ^1^) and 
in fig. [5] for the relaxation regime (i 3> i*), respectively. 
Given a specific force ratio, the evolution of AR\\ (t) corre- 
sponds to a vertical line through fig. [5] until i — t where 
it crosses over to the relaxation regimes {t ^ t^,) depicted 

in fig. [HI since in most cases ~ t^. 

Let us first turn to the propagation regimes of fig. [5] 
In the universal initial regime (light shaded) the scaling is 
linear in the force difference Af = /ext — /pre, see eq. p^ : 
it is followed by a quasi-static regime (white) with different 
force scaling for asymptotically small (<) and large (>) 
force ratio; in these limits, we find analytical values for 
the prcfactors b'' and as shown in the figure caption. 



B. Obermayer et al.: Stretching dynamics of semiflexible polymers 



9 




r{15/8)^P '^J^ 



/pre/Zext 

Fig. 5. Regimes of intermediate asymptotics (separated by 
thick black lines) for the bulk contribution to the change in 
end-to-end distance /i-Ry {t) (boxed formulas) in the propaga- 
tion regime f ^ time t/f~^^ vs. force ratio /pro//cxt (log- 
log scale). Asymptotic values for the prefactors are &^ ~ 2, 
6^ ~ —2.48, and a = 4(2/7r)^/*/-\/3. See text for explanation. 



(/pro = 0), the other with large force ratio connects to the 
"release" -case (/oxt = 0) of ref. [27], respectively. In the 
final regime t ^ with /ext,/pre ^ /c, the to leading 
order time-independent correction SF to the equilibrium 
tension profile [eq. (|26p ] implies via eq. ([5^ that ^R\\ (t) 
prematurely attains its equilibrium value 



.-1/2 
Jn 



I pre 



„-l/2 
J cxt 



(34) 



From eq. (|26f|) we find the subdominant correction term 

AfL 



(35) 



The intermediate regime of homogeneous tension re- 
laxationtj! ^ t ^ is discussed in ref. [27] and Ap- 
pendix lX^ For /pic 3> /c, /cxt we obtain from eqs. (|56bl32p 
the expansion: 



AR\i{t) = - 



1/3 



LO 



2/3 



(36) 




1 /pre/ /ext 

Fig. 6. Regimes of intermediate asymptotics for AR\i (t) as in 
fig- [SI but here for the relaxation regime t ^ ti,. See text for 
explanation. 

The general formula follows from eq. ([221) and reads 

AR\\ it) = ~idMOK'^'fl!tt'^', (33) 

where the dependence of the slope d^(p{0) on the force 
ratio c = /pre//ext IS depicted in fig. H) The intermedi- 
ate taut-string region (dark shaded) emerges only for very 
small force ratio. Here the analytical value for the prefac- 
tor a applies to the case /pre = 0. 

Except for the case /pre 3> /ext,/c, the propagation- 
relaxation crossover at t ~ i| leads directly to the regimes 
depicted in fig. [51 Here we find again a short-time regime 
linear in Af, followed by two different intermediate regimes: 
the one with small /pre//cxt connects to the "pulling" -case 



5.2 Numerical results 

Using the algorithm outlined in sec. [H we solved eq. ([291 
assuming a polymer of finite length L for various values 
of the two parameters c and A^. The observable Z\i?j[(t) is 
then extracted from the slope of the profile at s = (3 and 
s ~ L via eq. ([5^ . As we analyzed the effect of the ratio 
/prc//cxt of the two force scales present in our scenario in 
the previous sec. 15. 1|, we concentrate now on the magni- 
tude of the force /ext relative to the critical force /c via the 
control parameter = a/ /c//oxt for fixed ratio /pre//ext- 
In fig. [71 we depict results for /pre//ext = 10~^ and vari- 
ous values of A^, i.e., for a sudden increase in stretching 
force. One can discern the intermediate regime at times 
/cxt ^ t ^ (/ext /pre) which vaiiishcs if the stretch- 
ing force is changed only slightly. In the limit /pre = 0, it 

connects to the "nonlinear MSPT" propagation (t <^ i|^) 
and "nonlinear OPT" relaxation {t ^ t^) regimes of the 
"puffing" -scenario discussed in ref. [19], respectively. Fig. [51 
shows results for the inverse scenario, a sudden decrease in 
the stretching force, with /pre//cxt = 10''. Here, the time 
window /j^c ^ t <$:i produces a separate asymptotic 
regime only in the relaxation phase which connects to the 
regime of nonlinear OPT relaxation [19] of the "release" - 

scenario if /cxt = 0. The time window t| ^ t ^ for the 
regime of homogeneous tension relaxation with its t^/^- 
scaling gets broader as /pre //ext ^ co- 

Because of purely numerical noise in figs. [3 and [HI 
we have replaced the long-time asymptote of AR^^{t) for 
yl^ ^ 1 with dotted straight lines, because the subdomi- 
nant C'(<^^/^)-corrections in eq. ([SS]) become smaller than 
the numerical resolution and it is not possible anymore to 
reliably evaluate eq. ([5^ . The dashed lines in figs. [3 and 
m indicate the asymptotic growth laws for (t) in the 
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Fig. 7. Numerical results for AR\\{t) for /prc//cxt = 10~* and 
= io~* - * from top to bottom. The dark vertical lines mark 
the crossovers at t = f~Ji (left) and t — (/ext/pro)~^ (right), 
while the dark diagonal line denotes the crossover ai t — t^j^ ~ 
t*, see table [1] The dashed lines indicate the asymptotic power 
laws of figs. [5] and [6] 
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Fig. 8. Numerical results for AR\\(€) for /pre//ext = 10* and 
yl^ = io~* "* from top to bottom. The dark vertical lines mark 
the crossover at /; = /p^^ (left) and t — (right), the dark 
diagonals mark the crossover at t — t^^^ (upper) and t — t^, 
(lower) . See table [1] The dashed lines indicate the asymptotic 
power laws of figs. [5] and (6] and of eq. (|36|) . 



respective regimes summarized in figs. [5] and [51 By actu- 
ally collapsing the corresponding tension profiles onto the 
scaling functions derived in sec. [3] and appendix I A. 2) we 
find excellent quantitative confirmation of our analytical 
approximations in the suitable asymptotic limits (results 
not shown). 



5.3 Experimental implications 

The numerical results for AR\\{t) as shown in figs. [7] and 
[8] over 20 time decades confirm the intermediate asymp- 



totics of eq. (jTO|) . They can also be used to fit actual 
experiments, if one keeps in mind that these unbounded 
asymptotic growth laws arise from superponing the ex- 
ponential relaxation of (infinitely) many modes of wave- 
length < g < cxD with relaxation times Tq ~ q'"^ / {q"^ + f) , 
see the response function x±(9; t') [eq- (HI)]- A real poly- 
mer provides mode cutoffs at (/max — """/a where a is some 
microscopic length scale (say, the polymer thickness) and 
at (/mill — tt/L. The latter cutoff defines a time scale tj^, 
which is the relaxation time of the longest bending mode: 



ti 



1 

cxt 1 



if /cxt < L-^, 
if /cxt > L-^- 



(37) 



For times beyond tj^, the continuum limit L — *■ 00, i.e., 
disregarding the proper boundary condition and the dis- 
creteness of the mode spectrum, becomes invalid; the al- 
gebraic relaxation ends and is followed by normal expo- 
nential relaxation. Since the crossover times i*, and tj^ 
depend quite strongly on experimental control parameters 
such as stretching force, filament length, and persistence 
length, it is nevertheless possible to cover a broad win- 
dow of stretching and relaxation dynamics {t ^ tj^) by 
adjusting these parameters. Our numerical solutions can 
easily be adapted to the specific experimental setup at 
hand and provide a quantitative description without free 
parameters. 



6 Conclusions 

We have analyzed the generic scenario of a semiflexible 
filament prestretched with a force /pre that is suddenly 
changed to /ext- Related setups involving only one force 
scale have been discussed previously [13-16], but with di- 
verse results based on partially contradicting approxima- 
tions. Extending a systematic theory in terms of backbone 
tension [19] to include two force scales allows to resolve 
these discrepancies. Based on this tension formalism (for 
its detailed formulation see ref. [25]), we derived equations 
of motion and motivated an intuitive blob-picture in or- 
der to estimate growth laws for a central quantity, the 
boundary layer size or tension propagation length ^||(t), 
which is a measure for how far (or how fast) longitudinal 
correlations spread along the filament. The scaling laws 
[eq. p3)) ] for (t) have been confirmed by computing ten- 
sion profiles as scaling solutions to asymptotic differential 
equations. We find that, for times beyond the initial linear 
regime, which is dominated by the relaxation of very stiff 
bending modes, the prestretching force /pre has a signifi- 
cant influence. If /pre = 0, slowly relaxing long-wavelength 
contributions in the initial conformation resist equilibra- 
tion and dominate the relaxation process: the taut-string 
approximation applies. If /pre is finite, these modes are re- 
moved from the initial configuration. Hence, the dynam- 
ics equilibrates locally, which allows for a quasi-static ap- 
proximation. Only if /pre /ext is very small, there is 
a crossover from the former to the latter situation. The 
limiting single-force cases /cxt = and /pre = 0, respec- 
tively, analyzed in detail in ref. [27], were also recovered. 
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However, the associated hmits were found to be highly 
nontrivial, which is reflected quite subtly in the corre- 
sponding equations of motion for the tension. These equa- 
tions follow systematically from the general eq. (fTO|) . while 
the blob picture provides a more intuitive understanding 
of the underlying process and the complicated interme- 
diate asymptotics it produces. We presented a numerical 
algorithm to solve eq. (jlOp for very general initial and 
boundary conditions, derived new predictions for exper- 
imentally relevant observables and checked them against 
numerical solutions. Finally, the applicability of our re- 
sults to experiments was discussed. Our findings for the 
scenario analyzed in this work should be relevant for the 
viscoelastic response of more complex structures such as 
biopolymer networks [2,5]. Further, we argue that the 
above established dependence of the relaxation process on 
the initial conditions also generalizes to other force pro- 
tocols, such as (weakly) time-dependent external forces, 
transverse forces, elongational flows, or scenarios involv- 
ing sudden temperature changes. 
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A Analytical approach (details) 

A.l Asymptotes of eq. (fTOjl 

Following the analysis of ref. [27] , we note that the Fourier 
integrals in eq. (jlOp are dominated by wavenumbers near a 
(jni where the exponents are of order 1 . Explicitly, we find 
for given time t and integrated tension F the asymptotes: 

iiFyt<t:i (38a) 

^"^ ~ \f-^^^, if i^V* > 1 (38b) 

Observing that — scales like the inverse of the 
equilibration length £_l [eq. (fT2|) ]. we conclude that the 
first asymptote comprises free (bending-dominated) relax- 
ation where in the response function [eq. (jj])] the tension 
contribution q^F is subdominant, whereas the second one 
corresponds to forced (tension-driven) relaxation with the 
subdominant bending contribution q^t. 
Secondly, we write eq. pH)) as 



diF = D- S. 



(39) 



The two contributions creating curvature in the tension 
profile have a direct physical interpretation. 



D = 



""^P 



+ ./pro 



(40a) 



is of deterministic origin since it expresses the relaxation 
of initially excited modes for a polymer equilibrated under 



the force /pre. The second term 

r^2q'- /*dt'e-29'[«'(*-*')+^^(«.*)-^'(^^*')] (40b) 
Jo ""^P ^0 

is of stochastic origin since it comprises the accumulated 
influence of thermal excitations from times t' < t mediated 
with the response function x± [eq. (SI)]. With the general 
distinction [eq. (j38p ] in mind, we first turn our attention 
to the term D. In the bending-dominated case [eq. (|38ap ] 
with F^/i ^ 1, we neglect the tension contribution q^F 
in the exponent and obtain 

which possesses, depending on the magnitude of /pre, the 
two asymptotes 

93/4^1/4 

Dl.a ~ TTTTTTTT- if /pre* « 1, ^ « 1 (42a) 



F(l/4)^p 
1 

^^-pjpre 



if/pV»l, FVt«l. (42b) 



In the tension-driven regime F^/t S> 1, on the other hand, 
we get for D by neglecting the bending contribution q^t 



dq 1 - c-^""^ 

""^p + /pro 



ifFV^>l, (43) 



which can be expanded for asymptotically small or large 
ratios of the product Ffp^c- The leading-order terms are 



2F 



D2m~J^ ifi^/prc«l, F7i»l (44a) 



1 



2,6 



^^-p^prc 



ifi^/pro»l, F'A»1- (44b) 



The stochastic term S is simpler. In the regime F^/t ^ 
1 we neglect the tension terms and simplify 



Jo Jo 

dq l-e-2«'* 



(45) 



23/4^1/4 



if F^/t < 1. 



F(l/4)^p' 

For F'^/t ^ 1, the t'-integral in the stochastic term S is 
dominated by contributions near the upper limit and we 
can linearize F{s,t) - F{s,t') « dtF{s,t){t - t'): 

S^S2^ r^2q' rdt'e-2«'(*-*')[«'+*^] 
Jo ""^P Jo 



dg 1 — e 



"^P 
°° dq 1 
nip g2 + dtF 



dtF 



(46) 







2LV^' 



if F^/t > 1. 
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In the third line wc neglected the exponential in the limit 
J^V* > 1> which expresses the instantaneous equilibration 
underlying the quasi-static approximation of ref. [16]; this 
is an expansion similar to the one leading to cq. (j44bp . 



A. 2 Aymptotic equations of motion and solutions for 
the relaxation regime 

The asymptotic forms of eq. pO|) used in sec. [3] were mo- 
tivated physically and based on approximations employed 
previously [13,15,16], but can also systematically be de- 
rived. The right hand side D—S of cq. (|39p was analyzed in 
app. lA.lj for asymptotes of the deterministic and stochas- 
tic term D and S, respectively, see eqs. ()42I44I45I46|) . Now 
we can treat these various limits separately and find so- 
lutions for the relaxation regimes omitted in sec. [3l For 
t 3> we can exploit the fact that in this time regime 
the tension deviates only by a small correction 5F from 
the flat equilibrium profile /oxt^- It is therefore possible to 
expand the right hand side of cq. (fTOll about this constant, 
which to lowest order gives &l5F w [D — ^]F^f tf 

1. If F^/t < 1 and f^^J, < 1, we obtain D k. Di^a and 
S Si, which cancel to leading order. But since the 
preconditions imply small F and /pre, the right-hand- 
side of cq. ((39)) can be linearized with respect to these 
quantities. This was performed in sec. 13.11 see eq. (fT4|) . 
and lead via a Laplace transform to the tension profile 
eq. pT|) . In the relaxation regime L \ with \[z) = 
2^/^£y^z~^/^, we obtain the following tension profile: 



F[s,z) 



/o: 



s{s ~ L) 
z2 2A2 ■ 



(47) 



Backtransforming, we get 



F{s,t)^ U,t+-^j-^——s{s-L), (48) 



2V4r(7/4)^p' 
i.e., a small correction 5F to the equilibrium tension 

2. F'^/t < 1 and f^^^t > 1 gives D w L>i.b and S w Si, 
hence 



dtF = Di,, - Si 



-Si. 



(49) 



This implies constant curvature (no propagation) in 
the time regime <C i ^ ■ The tension reads 

9-1/4+1/4 _ T\ 

3. FVt > 1 and Ffp,^ < 1 yields 

~ D2,a -S2^ D2,a, (51) 



which comprises the taut-string case [eq. (|21a[) ]. We 
discussed the propagation regime in sec. (|3.2|1 : the re- 
laxation regime follows by writing F{s, t) = /oxt^ + 
SF{s,t) and expanding, which yields 



(52) 



4. F'^/t > 1 and Ffp^c > 1 leads to 

d^sF ~ D2.h - S2, (53) 
which corresponds to the quasi-static approximation 
[eq. (|23p ] . We discussed the propagation regime i <C t 
and the regime of premature tension saturation t ^ 
in sec. 13.31 In the latter case, we used for eq. (|26[p the 
next-to-leading order terms of the asymptotic expan- 
sions eqs. (|44bl46p . For /p^ > /ext and /pre » /c re- 
mains an intermediate regime t| ^ i <C of homoge- 
neous tension relaxation, see table [TJ It turns out to be 
very related to the "release" -case /ext = of ref. [27]. 
It is not possible to expand F{s, t) = /oxt t + SF, since 
the correction 6F is for t ^ ti, not (yet) small com- 
pared to /ext t. However, as in ref. [27] for /ext = 0, we 
can solve eq. ([M)) using the separation ansatz /(s, t) = 
g{t)h{s / L) and compute correction terms in the limit 
/pre 3> /ext- Disregarding the proper initial conditions 
and assuming Dirichlet boundary conditions for the 
moment, we find that asymptotically 



• r2\ 2/3 



(54) 



Although this is a decreasing function of time, wc find 
that if i ^ /(-p the preconditions {pt 1 and 
Ff pre ^1) that lead to the asymptotic cq. are 
fulfilled. In order to correct the Dirichlet boundary 
conditions assumed above to h{0) = h{l) = ^{t) = 

/oxt/5(t), we require that t ^ Lr' I {(-pftit ) which makes 
7 <C 1. Both conditions reduce to t ^ t*, cf. tabled] 
The roughly parabolic profile h{a) is a solution of 



h" 



1 



with h{Q) = h{l) = 7(i), (55) 



6/1I/2' 

where a subdominant time-dependence stems from the 
small correction term 7(t) to the Dirichlet boundary 
conditions. It has the following characteristics: 

Ml/2) = 3^(0'^' + i7 + O(7'/'), (56a) 

h'{Q) = 12-1/3 _ (2)2/3 ^1,2 ^ 0(^3/2)^ (5g^^) 

B Numerical approach (details) 

B.l Evaluation of the nonlinear term 

The computation of the nonlinear term G with its de- 
pendence on <P{a,Tk) requires to evaluate the z-integral 
in eq. (|29p . This can be done analytically, if a piecewise 
linear interpolation is used in the exponent: 



<Z>(cr,T„) - z<P{a, ZTn) ^ Ak + Bk{zk+i - z), 



(57) 



for an index fc G [0, n — 1] such that Tk < t„z < t^+i and 
with Zfe = TkjTn- The linearization coefficients are 

Au = <P{a,Tn) - Zk+l^i(7,Tk+l) 

<P{a,Tk+i)-^<J,Tk) (58) 



Bk =^ia,Tk+i) + Zk+i- 



Zk+l - Zk 
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With the linearization ([57|) . G'[^((t, t„)] is expressed as: 



fc=0 



1 _ Q-2fr„[q^+-P{cr,T„)] 



-2lJ^T„(Zfc + i-2fc)[lJ^+Bfc] _ 2^ 



(59) 



It turns out that the nonlinear problems [eq. (j3ip ] are 
ill-conditioned in many cases. Their solutions are strongly 
affected by numerical errors and often the algorithms used 
to solve these problems fail to converge. Hence, it is neces- 
sary to evaluate the nonlinear term [eq. (|59[) ] as accurately 
as possible. In order to represent the intermediate asymp- 
totic regimes, the variables r and <P have to vary over many 
orders of magnitude. Thus, in eq. we have to take 

care of overflow/underflow artefacts from the exponential 
functions. For the numerical evaluation of the g-integrals, 
we need to make sure that the dominant contribution from 
modes near (/m = gmV/oxt is properly taken into account; 
see eq. ((38)) . 



B.2 Nonlinear boundary value problems 

When solving a two-point boundary value problem with 
general mixed linear boundary conditions such as eq. (|3ip 
(where we omit the dependence on t„), the problem is 
mapped to a system of nonlinear equations. Writing the 
values ^((Trn) = ^m, whcrc m = 0, . . . , M, as vector 4* ~ 
(<?0j ■ ■ ■ ,^m)^ , the problem can be posed as 



T(#) = 0. 



(60) 



The nonlinear operator T is defined via the residuals of 
eq. dSIl): 



(T(#)),„ = A„*-G[<?,„], m = l,. 
(r(#))o = aoo^o + aoiD^^ - am 
{T{^))m = aio^M + aiiDl;^ - ai2 



,M -1 



(61) 



Here, we have used discrete representations of the differ- 
ential operators So- and 9^. For the possibly non- uniformly 
discretized coordinate (t,„ with local stepsize = cr,„+i — 
CFm and /i~ = am - cr,„_i, they are given by: 



2 



I I'm' I'm 



T+(<^1 -<^0) 



(62) 



In the formulation [cq. ([60]) ] our task is to find the 
zero of T in the vector space of discrete function rep- 
resentations The common strategies are all based on 
Newton's method [28]. Given an initial guess ^o, this it- 
erative method seeks a correction d such that + ^ is 
the solution. If the initial guess is good, the correction is 
small and wc can write 

r(#o + d) « T{^o) + J{^o)d = (63) 
with the Jacobian matrix 



(64) 



This matrix can be computed either analytically or by 
finite-difference approximation. In principle, the linear equa- 
tion [cq. (|63|l ] can be solved straightforwardly for the op- 
timal correction d: 



d = -J-i(<fo)T(<fo). 
Finally, by iterating 

=#fe- J-i(<ffc)T(«ffc), 



(65) 



(66) 



one should finally find the solution ^* = linifc^joo with 
r(^*) = 0. For "wcU-behavcd" operators T, this iteration 
converges quadratically in the vicinity of the solution. Un- 
fortunately, in our case T is often ill-conditioned near the 
solution. More sophisticated methods therefore modify the 
iteration [eq. ([5^ ] [28]. Still, a very accurate evaluation 
of the nonlinear term Gl'P] in eq. ([61]) is crucial to ensure 
convergence. For our implementation, two different algo- 
rithms in the GNU scientific librarjo and in the NLEQ 
package!! have been used. 



B.3 Error control 

Wc can control discretization errors stemming from rep- 
resenting the continuous problem [eq. ([29]) ] as system of 
nonlinear equations [eq. (|60p ] by reliable error estimates. 

An estimate for the spatial discretization error can 
be obtained via a deferred corrections approach [31]. The 
nonlinear system [eq. (|60p ] is written with the low-order 
accurate representation [cq. ((Ml)] for the spatial deriva- 
tives. Given its low-order accurate solution # and another, 
high-order accurate operator T*, we can solve the per- 
turbed system T{^*) — —T*{^) to obtain a high-order 
accurate solution . The discretization error of ^ is es- 
timated by the difference ^* — The important point 
is that not the (possibly complicated) high-order system 
associated with T* , but only the perturbed low-order sys- 
tem needs to be solved. A high-order accurate operator T* 
is obtained by using more accurate representations of the 



The GSL is available from www . gnu. org/ sof tware/gsl and 
uses Powell's hybrid algorithm [29]. 

^ The NLEQ package is available from www . zib . de and doc- 
umented in ref. [30]. 
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derivatives than those in eq. (|62p , in our case by means of 
Noumerov's method. Altogether, we obtain at least three 
digits of accuracy for our choice of discretization. We use 
"graded" grids, ct™ = aM{rn/MY for to = 0, . . . , M - 1. 
For the semiinfinite polymers, we take (3 = 1.25 and dy- 
namically adpat <tm] otherwise we take (3 = 1 and um ~ 1- 
The discretization error stemming from the r-discre- 
tization needs to be estimated because we want to cover 
many orders of magnitude in the time variable. Our re- 
sults show that an exponential growth in t„ does not in- 
crease the discretization error (this is due to the right hand 
side of eq. ()10p becoming asymptotically independent of 
shorttime effects in the nonlinear regime, cf. eq. (j46ll ). 
Hence, we choose t„ ~ Z\((5" — l)/((5 — 1) for n = 0, . . . , iV 
with the initial stepsize A and the growth factor 5 > 1. 
To estimate the discretization error, we take "fine" steps 
•f„ = ro((5"/2 - 1)/{S - 1) with n ^Q,...,2N such that 
T2n = T„- In the step n two additional "fine" solutions 
'^i<7,T2n-i) and ^((7, f2n) are computed. The error esti- 
mate follows from the difference ^(cr, f2„) — <?(cr, t„). Al- 
together, we get at least three digits of accuracy for our 
choice of the crucial parameter S: for the semiinfinite poly- 
mers, we take S = 1.1 (only tension propagation), other- 
wise d = 1.05 (also tension relaxation). 
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